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Numerical relativity codes now being developed will evolve initial data representing colliding 
black holes at a relatively late stage in the collision. The choice of initial data used for code devel- 
opment has been made on the basis of mathematical definitiveness and usefulness for computational 
implementation. By using the "particle limit" (the limit of an extreme ratio of masses of colliding 
holes) we recently showed that the standard choice is not a good representation of astrophysically 
generated initial data. Here we show that, for the particle limit, there is a very simple alterna- 
tive choice that appears to give excellent results. That choice, "convective" initial data is, roughly 
speaking, equivalent to the start of a time sequence of parameterized solutions of the Hamiltonian 
constraint; for a particle in circular orbit, it is the initial data of the steady state solution on any 
hypersurface. The implementation of related schemes for equal mass holes is discussed. 



I. INTRODUCTION 



Recently the numerical computation of gravitational waves from coalescing black holes has received considerable 
attention and effort for closely related reasons, (i) The problem is a testbed for numerical relativity and has been the 
center of a collaboration of eight universities [Q. (ii) The results are necessary to understand details, and perhaps even 
gross features of coalescing black holes, presently viewed by many as the most promising source of gravitational waves 
detectable by earth based instruments ||,^. (iii) The problem is interesting in its own regard as the first solution of 
Einstein's equations for strong field interactions without simplifying symmetry. 

The job of the numerical codes under development will be to evolve forward in time from an initial value solution 
on a starting hypersurface. There are many ways that initial data can be chosen to represent the starting state of two 
black holes. If the black holes are far apart on the initial hypersurface, then any reasonable choice of hypersurface 
data will suffice. Unfortunately the instabilities in evolution codes mean that the codes can only run for a relatively 
short time, and initial data will have to be specified on a hypersurface as late as possible. 

For the most part, the initial data used to the present has been that put forth by York and Bowen, and coworkers 
[Q-^. This scheme starts by requiring that the initial geometry be conformally flat and that the extrinsic curvature 
have a property called "longitudinal," equivalent to the condition that it be derivable from a vector potential. The 
solution is then made inversion symmetric, so that topological throats representing the black holes connect two 
isometric asymptotically flat universes. (This step of making the hypersurface data inversion symmetric turns out to 
destroy "longitudinality" [Q. ) 

The advantages of this approach to initial data include the following: (i) Since it is widely adopted in the numerical 
relativity community, comparison of results of different groups is facilitated, (ii) It is specific and definitive and the 
numerical methods of solving for these initial data are highly developed (iii) The inversion symmetry gave a 
simple inner boundary condition for evolution codes, (iv) An argument could be made that the initial data was free 
of, or had a minimal content of, radiation encoded (in some sense) in the initial data themselves, so that the outgoing 
radiation found in a spacetime evolution should represent only the radiation due to the dynamical interaction of the 
black holes. 

The Bowen- York (hereafter BY) scheme served very well in the role intended for it. One should ask, however, 
whether it is the appropriate scheme for computations supporting the search for gravitational waves. There are 
strong reasons for concluding that the answer is "no." One immediate reason is that the Kerr geometry itself cannot 
be put in this scheme; a constant time slice of the Kerr geometry is not conformally flat. Even a single spinning BY 
hole will therefore radiate gravitational waves as it settles down into its Kerr final state ^ . 

A second rather direct reason comes from recent calculations we carried out for radiation from a particle falling 
radially into a Schwarzschild hole ||l^,|ll|, hereafter Paper I and Paper II. In these calculations, the "particle" was 
taken as a perturbation of the background black hole spacetime and computations could be done relatively easily 
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for a variety of choices of data on an initial hypersurface. Furthermore, evolved data could be extracted on later 
hypersurfaces and compared with various prescriptions (BY etc.) for hypersurface data. In the approach used, based 
on the perturbation formalism of Moncrief , the initial and evolved data was represented by a spacetime function 
ip that is totally gauge invariant, so that the interpretation of the physics was not tainted by questions of coordinate 
choices. 

The calculations were used to study a set of hypersurface data that somewhat generalize the BY data. Those 
sets considered were all conformally flat and longitudinal (CFL). In practice this gave, for each multipole of the 
perturbation, a complete set of differential equations that had to be satisfied by the hypersurface geometry and 
extrinsic curvature (or, in practice, ip and its time derivative ip). The only remaining freedom is the choice of 
boundary conditions (consistent with asymptotic flatness). These choices are equivalent to the choice of conditions 
deep inside the throats representing the holes, i.e., in the asymptotically flat universes at the end of the throats away 
from "our" asymptotically flat universe. We focused on two choices of these boundary conditions. Once choice was 
the choice of the BY data itself. The second was the choice that the solution at the horizon be "frozen" or "matched" 
to its initial value. This approach was based on the fact that very near the background horizon a simple analytic 
solution could be given for the evolution of the data on the initial hypersurface to subsequent hypersurfaces. The 
solution for ip was simply that it was "frozen" to retain its initial value at the horizon. The solution for ip reflected 
more detail, and was said to be "matched" to the initial solution at the horizon. We could also choose to have ip 
"frozen," but we could choose the boundary conditions for ip to be the BY conditions, and so forth. 

For a particle falling from rest at a large initial distance, we found that the solution that evolved, and the CFL data 
on later hypersurfaces were in serious disagreement. If evolved data were replaced by CFL data on a late hypersurface 
the resulting waveforms and radiated energies were drastically altered. We found furthermore: (i) In the class of 
CFL hypersurface data the best choice is the "frozen-matched" choice, but that the difference between this and other 
CFL choices is small compared to the disagreement with the evolved solution, (ii) The disagreement between evolved 
hypersurface data and CFL data does not lie, in any obvious way, in the "radiation content" of the data, but rather 
in the details of the data near the location of the particle, (iii) The disagreement between the true and CFL data was 
much more severe in ip than in ip. 

In the present paper we look again at the problem of head on collision of two nonspinning holes in the particle limit, 
and ask whether there is a more successful alternative prescription for hypersurface data. More specifically, let us 
consider a configuration of black holes at late times that evolved from a precursor configuration at large separation, 
and let us suppose that we know the positions and momenta of the holes on the late time hypersurface. Is there a 
prescription for the data on the late hypersurface which is in good agreement with the "true" hypersurface data (i.e., 
the data evolved from the large-separation precursor) , and which gives an accurate prediction of the radiated energy 
and waveform? The answer to this question seems to be that there is, and the prescription, at least in the particle 
limit, is a very simple one: that the data be "convective," in a sense to be defined below. 

In Sec. IIA we sketch the mathematics of the Moncrief formalism we use and of our notation and conventions, and 
we present some results of evolved data vs. CFL data; this sketch will be quite brief; a reader interested in the details 
will find them in Paper II. The proposal for convective data is discussed in general in Sec. IIB, and is implemented in 
a particular manner, for radial infall, in Sec. IIC. Numerical results from this implementation are given in Sec. IID. 
A discussion of the implications and possible extensions of this work is presented in Sec. III. 



II. ASTROPHYSICAL INITIAL DATA 



A. Review of particle results for BY initial data 



1. Moncrief-Zerilh formalism 



The numerical results we will present below will be for a particle infalling radially into a Schwarzschild black hole. 
Due to the axial symmetry of our problem, and the absence of rotation, the perturbations are pure even parity. 
To describe the perturbations, we use the Regge- Wheeler |Q notation, but not the Regge- Wheeler gauge. In this 
notation the even parity metric perturbations are decomposed into scalar, vector and tensor spherical harmonics, 
with coefficient functions Hq, Hf, hQ,h{, K^, , depending on r, t the standard radial and time coordinates of the 
Schwarzschild background. For simplicity, we drop the multipole index £ on these perturbation functions. 

The perturbation coefficients H2,hi, K,G describe only the perturbed 3-geometry, and are independent of the 
choice of shift and lapse. In terms of these, Moncrief ||l2| defines the totally gauge invariant function 



ip{r,t) 



A + 1 



^ + - rdK/dr} 

\r + 6M 



+ -(r - 2M) {r^dG/dr ~ 2hi) 



(2.1) 



2 



where we have used Zerilh's normaHzation for tp and his notation 

X={e + 2)(£-l)/2. (2.2) 
The basic wave equation for an infaUing particle is given in Paper I as 

Here r* = r + 2A/ln(r/2Af — 1) is the Regge- Wheeler "tortoise" coordinate and Vi is the Zerilli potential 

^^('^^ rHXr + 3My " ^'''^ 

The right hand side of Eq. (2.3) is a source term constructed from the particle stress energy 0Jll[. For a point 
particle of proper mass mg, the stress energy is given by 

T^"" = {mo/U")U^'U''6[r - rp{t)]S^[n]/r'^ , (2.5) 

where is the particle 4- velocity, and rp{t) gives the radial position of the particle as a function of coordinate time, 
starting from rp ~ vq at t = 0. 

The total energy radiated in the ^th multipole, after a hypersurface at to is 



2. Functional freedom in the hypersurface data 

An even parity multipole perturbation of the 3-geometry has the four functional (of r, t) degrees of freedom contained 
in hi, K,G, but these must satisfy the (perturbed, multipole) Hamiltonian constraint, a single second order 
differential equation containing the particle source. The choice of a conformally flat 3-geometry reduces the problem to 
a single functional degree of freedom (since, in an appropriate coordinate system, hi = G = Q and H2 = K) satisfying 
the single Hamiltonian constraint. The only freedom is the choice of the constants specifying the homogeneous solution 
of the Hamiltonian constraint. One of those constants is fixed by the condition that the solution be well behaved 
at spatial infinity. The remaining constant can be specified to "freeze" the horizon value of ■(/' at its initial value, to 
make the solution inversion symmetric, or in other ways. If the requirement of conformal flatness is dropped, then 
any sufficiently smooth function can be added to a particular solution ij: that solves the Hamiltonian constraint. 

For -0 things are a bit more complicated. For an even parity perturbative extrinsic curvature there are four functions 
of r,t, analogous to the four functions hi, H2, K, G specifying the 3-geometry perturbations. There are two second 
order equations, equivalent to Gtr = SnTtr and Gte = 0. If we require that the extrinsic curvature be longitudinal, 
then the extrinsic curvature must be derivable from a vector potential W, and for a given even parity multipole 
the functional degrees of freedom are W^{t,r) and W^{t,r). Since these functions must satisfy two second order 
differential equations, we have no remaining functional freedom. We have only to choose the four constants that 
specify how much our solution contains of the four homogeneous solutions of the momentum constraints. Two of 
these constants are fixed by the requirement that the extrinsic geometry be well behaved at spatial infinity. The 
remaining two constants can be chosen to make the solution of the BY type, to make the solution "match" the initial 
solution at the horizon, etc. If the requirement of longitudinality is dropped, then any sufficiently smooth function 
can be added to any particular solution of the momentum constraints. 



3. Results for CFL hypersurface data 



Here we briefly review the numerical results found for the CFL prescription. In the results presented here a particle 
started its infall from Schwarzschild radial coordinate Vp = = 15 (2M), where M is the mass of the spacetime (in 
the usual c = G = 1 units). At subsequent times, when the radial location of the particle is smaller, and the time 
label of a t = constant hypersurface is larger, we can extract the evolved ip, replace it with a prescribed CFL etc. 
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Figure |l] shows the result that is most important. At various times (parameterized by the particle location r* at 

these times) the numerical evolution was stopped, and the evolved quadrupole and ip were replaced by the CFL 
prescribed data representing the particle position and momentum on this hypersurface. Evolution forward from this 
hypersurface is then done numerically. The resulting energy is then numerically computed. The true £ = 2 energy 
for infall from tq = 15(2M) is 0.0164(m§/2M). If we evolve the prescribed data from a later surface we see that we 
get reasonably good agreement with the true energy for r* larger than around 5(2M), corresponding to Vp around 
4(2M). For smaller the use of prescribed CFL data badly overestimates the radiated energy. 

Comparisons of the true (i.e., evolved) quadrupolar tpy'ijj and the prescribed quadrupolar V': V' shown in the 
next two figures. Both are for a particle starting from rest at initial separation tq = 15(2M) at t = 0. For both, the 
comparison is made at t = 96.9(2M). This corresponds to the particle being at — 2(2M), and having momentum 
P = 1.27mo. A comparison of ^p is given in Fig. ^ in which the solid curve shows the evolved ip and the dashed curve 
is the CF tp for a particle of mass mg, at location Vp = 2(2A/), with momentum P = 1.27r7io. The corresponding 
comparison for ip is shown in Fig. ^. From these comparisons it appears that the CF solutions give a reasonably good 
approximation to ip near the particle, but in the strong field region of the background, the CFL form of V' has a very 
different shape near the particle than that of the evolved ip- We note here that the evolved ip will, in general, have a 
S[r — rp] behavior at the location of the particle. This behavior is exactly included in both the CFL and the convective 
(see below) prescriptions. 



B. Proposed astrophysical initial data 

Suppose that our particle source starts at some position ro,0O;Vo- On some later hypersurface, let us call it the 
t = ts hypersurface, suppose that we know the coordinate position xj^ and the momentum P^ for the particle. (Here 
latin superscripts, j, k, . . . are spatial indices running from 1 to 3.) From the geodesic motion of the particle we can 
find x^p and P^ at any time t. Let us now take ^{x^;Xp) to be a form of ip that is "appropriate" to the position (and 
momentum) of the particle. This could be, for example, the ip corresponding to a CF 3-geometry for the particle. The 
next step is to promote this solution to a function of time, by putting in the explicit time dependence Xp[t]. (Note 
that the resulting time-dependent function ^'(a;'^; Xp[t]) is not based on a solution of Einstein's equations, and will not 
solve the Zerilli equation.) Our proposal is to take the hypersurface data to be: 

iPU^=^{x'^;x';M) (2.7) 
<9*(x'=;a;^[t]) 



With this choice we are, in effect, saying that ip is changing in time only in order to adjust to the new particle position 
and momentum. The solution, in a sense, is dragged or transported to a new position. We call this the "convective" 
choice of ip. 

Initial data, of course, must be chosen so that it solves the constraint equations. But here (so far) we have been 
discussing the choice only of the single gauge invariant function ip on the hypersurface, and its time derivative. Since 
the constraints leave a functional degree of freedom in the initial 3-geometry and extrinsic curvature, ip and ip can be 
freely specified. 

One motivation for this choice of hypersurface data is that the disagreement between evolved and prescribed initial 
data seems to be rooted in the CFL ip much more than in ip. The proposed convective choice of hypersurface data 
allows us to maintain the CF form of ip, but find a new form of ip. A second motivation has to do with motion in 
which the hypersurface data is "obvious," a particle in circular orbit at frequency around a hole. In this case the 
appropriate approximation (neglecting radiation reaction to lowest order) is for the spacetime geometry to be periodic, 
with period 27r/il. More specifically, in terms of the coordinates t,r,9,(p of the Schwarzschild (or Kerr) background, 
all spacetime functions should have the form f{t, r, 9, (p) = F{r, 9,(p — fit). 

If we denote the position of the particle by rp — rg, 6p — 7r/2, ipp = ipo + i^t. The solution for ip must therefore have 
the functional form as 

iP{t, r, e, = $(r, 0, - nt) = $(r, 0, - ^p[t\ + (po) • (2.9) 

We can therefore take 

4'(x^a;'=[<]) = $(r,0,(^-(pp[t] + (^o) • (2.10) 



dxl[t] d^{x'';xl[t]) 



dt 



dxp 



(2.8) 
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This sohition clearly obeys Eq. ( p.SD on any hypersurface is, since the only time dependence in comes through its 
dependence on </'p[t]. If on some hypersurface i = is, one chooses initial data 



dip 



(2.11) 



then, in principle, numerical evolution would give a purely periodic solution. For x^^P^ corresponding to circular 
motion, one could, on the other hand, give other hypersurface data, such as CFL. The result would be a burst of 
initial radiation as the spacetime "settles down" to the periodic solution. If the circular orbit is meant to represent a 
late stage in a binary coalescence which is adiabatically decaying due to radiation reaction, then the periodic solution 
is obviously the appropriate one, and an initial burst of radiation is an anomaly due to the wrong choice of initial 
conditions. 



C. Convective conformally flat data 

We now apply the idea of convective ■)/' to a particle of mass mo falling radially into a Schwarzschild hole. As the 
appropriate form of V' on a hypersurface we choose a conformally flat solution, with "frozen" horizon conditions. That 
is, on any hypersurface, the r* — > —oo limit of ip is the same as on the i = hypersurface. 

The meaning of this solution can be nicely described with the notation developed in Paper II. For even parity £-pole 
perturbations of the 3-geometry, there are only two gauge invariant functions of r. One is the Moncrief ip defined 



by Eq. (2.1). The second is /conf, a gauge invariant measure of deviation of the 3-geometry from conformal flatness. 
(In the Regge- Wheeler [|l3| gauge /conf = H2 — K.) For all CFL choices /conf is zero on the hypersurface. This 
choice exhausts the functional freedom on the hypersurface, and the 3-geometry is then fixed by the Hamiltonian 
constraint. In practice this means that ^ satisfies a second-order differential equation in the r variable; the frozen 
horizon condition fixes the solution to that equation. 

We are then choosing the same hypersurface ■0 as a CFL choice. What is crucial is the difference in our choice of 
■0. Here our convective choice is tantamount to taking '^[[x^; x'^{t)\ to be conformally flat for all t. This corresponds 
to a form of V', on each hypersurface, with /conf = 0. It follows that 



^conf 



. (2.12) 



In order to understand better the meaning of this condition it is useful to relate it to the Hamiltonian constraint. To 
derive this relationship conveniently we momentarily specialize to the gauge /iq = /ii = G = of Regge and Wheeler 
In this gauge the Hamiltonian constraint is given by Zerilli ||lj] as his Eq. (C7a). For a particle at radial position 



Tp this takes the form 



2M\d^K / 5M\ldK A A + 2 / 2M\ldH^, 
1 „ o + 3 -— -K T^Ho — 1 — 



/ 2 A f \ 1 
[/o 1 „ ^ (2.13) 



where k = 2 mQ\/ A tt{2£ + 1), mo is the mass of the particle , an d [/° is its time component of 4- velocity. 



From Eq. (2.13) and its radial derivative, and from Eq. (2.1) we can write 



_ „ [(A - 3)r + 9M]M [27M^ + 24AMV + 3A(3A + l)Mr^ + 2A^(A + ly] , 

W-(r 2M)^ + (;^^^3^^)^ ^ (Ar + 3M)2r2 

nU'{l-2M/r)[X(X + iy-3M'] _ nU\r^2Mf _ 

+ (A + l)(Ar + 3A/)2 '^"^ "^^^ {X + l){Xr + M'iy^ ' ^ > 

This result was derived in the Regge- Wheeler gauge as an intermediate step, but is gauge- invariant. In a similar way, 
from the momentum constraints (see Eqs. (C7b), (C7d) of Zerilli fl^ ) and their r-derivatives, we obtain 



• , [(A-3)r-h9M]A/ ;, [27A/3 -|- 24AA/2r + 3A(3A + l)Mr2 + 2A2 (A -M)r3] • 
/conf = {r- 2A/)V. + (;,,^3^,,), [Xr + iM fr^ 

2k ?7°fp[9A/2 + 2A/(A - 3)r - A(A + 2,^] 



(A-|-l)r2(Ar + 3Af)2 



5[r - Tp] (2.15) 



KC/°rp(r-2A/)[9A/2-F2A/Ar-A(A + l)r2] KU^fp{r-2Mf 

(A + l)r(Ar + 3A/)2 ^'^ " ""^^ + (A -f l)(Ar + 3A/) " "^^^ 
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Note that this equation can be obtained as the time derivative of Eq. (2.14), but here we wish to emphasize its 
connection with the momentum constraints. The fact that there is consistency of the time derivative of the Hamiltonian 
constraint, and the s patia l derivatives of the momentum constraints, is simply an expression of the fact that Gq'' = 0. 

We can view Eq. ( 2.14 ) as telling us that /conf vanishes, since the right hand sid e van ishes for conformally flat data 
that satisfies the Hamiltonian constraint. Alternatively we can note that in Eq. (2.14) if we put /conf to zero, what 
results is a second order radial differential equation for Tp. In this sense, the condition /conf = and the Hamiltonian 
constraint gives us (aside from boundary conditions) the initial value of ip. In the same manner we can note that 
/conf = in Eq. (2.15) gives us a radial differential equation for ip. From this point of view the condition /conf = 
and the momentum constraints are fixing the form of ^p. By contrast, in the CFL prescription, the form of ip is fixed 
by the momentum constraints and the condition of longitudinality, roughly the condition that the conformal part of 
the extrinsic curvature can be derived from a vector potential. 

To see how our proposal works, let us first consider the radial infall of a particle into a Schwarzschild hole. We have 
seen in Paper II that there is a preferred choice, V'/: that keeps its value at the event horizon "frozen" 



2m(rp) V47r/(2^+l) ^ 
ipf = ^ ^ ^ \ „ I o a r — ^yr/r X 



A + 1 Ar + 3A/ 



(a + 1 + ^ + - 2M/r + ^) ^ r,„„(rp) ^^^p^ 

(a + 1 + M/r + ^1 - 2M/r [l + ^/¥JPj ) (fp/f)^ 
(^A + 1 + M/r + ^1 - 2M/r ( ^^l^vf^^ 



where 



f = (y^+^Jr-2M^ /\ , 
and where the upper expressions apply in the case r > rp, and the lower for f < fp. Here, as in Paper II 



(i + rw(r-o))-i 



(2.16) 



(2.17) 



(2.18) 



where Thomiro) represents the arbitrariness in choosing the data on the initial hypcrsurfacc t = (when the particle 
started at rest and was located at rp = tq), and 



m(rp) = imo 1 + 



1 



^1 - 2A//rp 



£0 



(2.19) 



where eq = -^/l — 2M/ro for a particle starting at rest when r — ro, and eo = 1/ ^1 — v?^ for a particle with velocity 
Voc at infinity, and C/° = eo/(l — 2M/rp). 

From this expressions we compute the "convective" ipc — T'pdrp4'f 



■ _ 2m{rp)^/I^/{2e+l) r.JVJf fp 



A + 1 



Ar + 3A/ [rp - 2A/) 



(^A + 1 + y + - 2M/r (l + ^/fj^^ ^ {{i + l)y^l -2A//rp + M/{rp + rp^\-2M/rp) 



(A//2) 



2^+1 



(2.20) 



(a + 1 + M/r + - 2M/r (l + v/^)) (l^/Y^2MJ7p - M/{rp + rp^l - 2A//rp)) (fp/f)^ 
A + 1 + M/r + ^1 - 2M/r ( " ^ - l)) ((^ + l)^/l-2M/rp + M/{rp + rp^l - 2A//rp)) {f/fpY+^ 



where the upper expressions apply in the case r > rp, and the lower for r < 



For radial infall f„ 



-(1 



2M/rp)^2M/rp + eg — 1/eo. Note that the choice of the freezing condition for ipf ensures that tpc vanishes on the 
horizon. This would not be true if we had taken either the particle limit of the Misner ||l^ or Brill-Lindquist p^ ] 
solutions for ?/;. It is also worth stressing that ipf^ does not depend on VhomiTo)- 

In principle, it could turn out that the CFL prescription and the convective prescription are the same. In Fig. ^ 
we show that this is not the case. This figure shows /conf for a particle dropped from rest at tq — 15(2A/). For the 
t = hypersurface, with the particle at rp — 15(2A/) (equivalent to r* — 17.64(2A/)), the solution is momentarily 

stationary, and hence /conf is momentarily stationary, that is, /conf = 0. Plots are given of /conf on hypersurfaces a 
short time after the initial hypersurface, constant-f hypersurfaces corresponding to particle positions r* — 14.55(2A/) 
and r* = 11.15(2A/). The figure shows clearly that /conf does not vanish. 
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D. Numerical results 



For some purposes the most important question to ask about a prescription for late time hypersurface data is 
whether it leads to a reasonably accurate estimate of the radiated gravitational wave energy that is generated in an 
astrophysical event. For the problem of a particle falling in from rest at tq, at time t = 0, we know with considerable 
accuracy the energy radiated in each multipole during the infall. We now must ask: if we replace the evolved data on 
a late hypersurface by some prescription (CFL, convective,. . . ), what radiated energy will we calculate? A numerical 
answer to this question is shown in Fig. |^. For a particle falling from rg = 15(2^/), the "true" quadrupole gravitational 
radiation emitted is 1.64 x IO^^itlq/ {2M). Figure || shows the radiation that is computed if the evolved data on a 
hypersurface is replaced with prescribed data and evolved forward in time. The prescribed data used for the 3- 
geometry is conformally flat with a frozen horizon boundary condition. The value of ip is given by the convective 
prescription. The results show that the use of such prescribed data, even at fairly late times, gives a very good 
approximation to the radiated energy. The worst case is an overestimate of the energy by a factor of around 2.4. This 
is to be contrasted with energy computed for prescribed CFL data, shown in Fig.^ where the error grows without 
bound as the hypersurface is taken to late times. The result is significantly in error only for data prescribed on a 
hypersurface when the particle radial position is around the peak of the Zerilli potential in Eq. (2.4), where we 
would expect the largest deviations in results based on any simple prescription. To clarify the relationship of the 
Zerilli potential and errors in the computed energy, the Zerilli potential is also indicated in Fig. |^. 

The success of the new prescription is demonstrated in Figs. ^ and ^ These figures compare the new prescription 
(conformally flat, frozen-horizon, convective) for hypersurface data with the true (i.e., evolved) data and with the 
BY data, on several hypersurfaces. In Figs. p^-R we use BY data as our example of the CFL prescription, since BY is 
the most familiar case. There would be negligible differences in appearance if we used instead any of the other CFL 
variations studied in Paper II. 

Figures ^ shows dip/dt for a particle falling from rest from position rg — 1.5(2Af) at time t = 0. This is not a 
reasonable astrophysical scenario, but it serves to emphasize strong field effects and magnify differences in hypersurface 
data. The improvement of the new prescription over BY data is dramatic, even on a hypersurface {t = 2.2[2M],rp = 
1.34[2M]) soon after the initial hypersurface. The new prescription agrees very well with the evolved data, both in its 
general shape near the particle, and in its numerical values; by contrast the BY data disagrees sharply with the evolved 
data. This disagreement grows at later times, and at t = 4.5(2M) the BY and evolved curves have little similarity. 
On the other hand, the new prescription gives a curve that agrees qualitatively with the true data, and agrees in size 
within a factor of two. The same comparison is made in Figs. for a particle falling from rest at ro = 15(2M) at 
t = 0. Here we see that both the new prescription and BY data agree with the evolved data reasonably well until the 
particle approaches the strong field region near the peak of the Zerilli potential. At smaller radii there are significant 
quantitative differences between the data of the new prescription and the evolved data, but these differences are much 
milder than the wide divergence of the BY data from the evolved data. 

A comparison of outgoing waveforms appears in Figs. |^. The solid curves in both (a) and (6) show the true waveform, 
seen by an observer at robs = 1000(2Af), for the quadrupolar part of ip, for a particle falling from rest, at time t = 
from position Vp = 15(2M) (equivalent to rp — 17.64[2M]). The dotted curve in (a) shows the outgoing waveform 
that results if conformally flat, frozen-horizon, convective data is substituted for true data on the t = 63.9(2A/) 
hypersurface. Similarly the dashed curve in (a) shows the waveform if the substitution is made at i = 93.5(2M). 
The three curves show remarkable agreement after the start of quasinormal ringing. Only at early times is there a 
noticeable difference: the long shallow dip in the true waveform is reduced or missing from the waveforms generated by 
the data of the new prescription. The dotted curve in (&) shows a waveform comparison at a yet later time. Although 
the new waveform has roughly the same magnitude of quasinormal ringing as the true waveform, it has a large phase 
shift relative to the true curve. This can be understood with the idea that quasinormal ringing is generated near 
the peak of the curvature potential. For the dotted curve in (&), the particle is placed at a smaller radius than the 
potential peak, and has no history of having moved inward through the potential peak. The quasinormal ringing 
excited by the prescribed data, then, must start at the time of the hypersurface on which prescribed data is defined. 
For the true curve, the quasinormal ringing starts earlier, at the time the particle passed going inward. Numerical 
experiments confirm with reasonable accuracy that the phase shift depends on particle position in this manner. 

In the work of Abrahams and Cook ||l7|, and Baker and Li the "close limit" was used to evolve prescribed data 
for the head on collision of equal mass holes. The close limit approximation applies to a hypersurface late enough that 
a single horizon surrounds both colliding holes, and only the large radius fields of the colliding holes are of importance. 
We have seen in Paper II, however, that in the particle limit this close limit does not work at any separation of the 
holes, if BY prescribed data is used. The method fails at large separation, because the close limit conditions fail; the 
method fails at small separation because the BY data is a bad representation of the true data. In Fig. ^ we see that 
the close limit does work at sufficiently small separations if the new prescription (conformally fiat, horizon frozen, 
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convective) is used. 

III. CONCLUSIONS 

We have seen that the convective prescription for initial data for a black hole coalescence gives a much better 
description of the way in which hypersurface data evolves than does the CFL prescription. The demonstration, 
however, has been confined to an example of a coalescence which (i) was a head on collision (ii) used a conformally 
flat 3-geometry, and (iii) involved the particle limit of nonrotating holes. 

It will be relatively simple to study another particle-limit example that is free of the first two of these limitations. 
As already mentioned in Sec. IIB, the steady state solution for a particle in circular orbit around a hole is an example 
of a solution that is convective for any time slice. It is straightforward to find this convective periodic solution in 
the form of a Fourier series, and this solution will not be conformally flat. This opens the possibility for several 
interesting comparisons. It will be reasonably straightforward to start with CFL data for a particle which at one 
instant has the correct position and momentum for a circular orbit. Similarly, it is simple to find initial conditions 
corresponding to a convective, initially conformally fiat data. Both these solutions could be evolved forward in time 
and the anomalous radiation could be found that is emitted as the field settles into its steady state solution. The 
amount of radiation energy in this anomalous initial burst is an indicator of how sensitive conclusions will be to the 
details of the initial data prescribed. If the energy is significant compared to the steady state emission during a single 
orbit, it is a sign that CFL initial data will give misleading results for black hole coalescence, at least for coalescing 
holes of very different mass. 

A considerably more difficult issue is the particle- limit context of our development. How are these ideas to be applied 
to coalescence of holes of comparable mass? In a very general way one might try to use the idea of a time sequence of 
3-geometries in order to fix the initial extrinsic curvature. There are many caveats that go with such a proposal, along 
with the question of just how one implements such a scheme? One could, for example, add time dependence to the 
conformal factor with the convective description, by making the parameters (location and momentum of the holes) 
time dependent. But we do not know at the outset what the correct dynamical sequence is for strong field motion. 
There is a more serious difiiculty. Such a scheme would give a definitive prescription for finding the time derivatives 
of the 3-geometries, and (with some physically motivated choice of lapse and shift) would give us the initial extrinsic 
curvature. This extrinsic curvature should in general give a result that does not satisfy the momentum constraint. 

A more promising approach might be to try to generalize the condition /conf = 0. One might, for example, use 
the Cotton- York tensor (3ij, which vanishes if and only if the 3-geometry is conformally fiat. By setting the 
time derivative of this tensor to zero one would capture, in the nonlinear theory, the same condition as /conf = 
in perturbation theory. A scheme like this, however, is based on conformal flatness, and this raises the question of 
whether conformal flatness itself should be retained in an astrophysically motivated specification of initial data. 

Somewhat related to the question of a convective prescription, but somewhat independent, is the question of the 
frozen horizon condition. For representing initial data that evolved from an earlier astrophysical configuration, the 
frozen horizon condition is "obviously" a correct constraint on data on a late hypersurface. Whether or not the 
convective prescription is used in connection with numerical relativity, the frozen horizon condition should be used if 
one wants a representation of an astrophysical problem. 
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FIG. 1. Radiated energy for replacement of evolved data by prescribed data. The specific prescribed data used was (CFL) 
conformally flat, longitudinal data with ^ at the horizon taken to be "frozen" to the starting value, and ip "matched" to the 
starting form. A particle was started from rest at ro = 15(2M) with momentarily stationary data, and the solution for ^ was 
evolved numerically. On various constant time hypersurfaces, labeled with the r* position of the particle on the hypersurface, 
the evolved data was replaced by the frozen-matched prescribed data, and evolved forward. The radiated i = 2 energy is 
plotted as a function of the hypersurface label r* at which the replacement was made. 

FIG. 2. Comparison of evolved and prescribed tp for infall from rest at ro = 15(2Af ) at t = 0. The comparison is shown for 
the hypersurface a,t t = 96.9(2M), at which time the particle has fallen to Vp — 2{2M) and has a momentum P = 1.27mo. The 
prescribed ip shown is for CF 3-geometry, with a frozen boundary condition at the horizon. 

FIG. 3. Comparison of evolved and prescribed ip for infall from rest at ro = 15(2M) at t = 0. The comparison is shown for 
the hypersurface at t = 96.9(2M), at which time the particle has fallen to rp = 2(2M) and has a momentum P = 1.27mo. The 
prescribed shown is for CFL extrinsic curvature with the boundary condition at the horizon "matched" to the conditions at 
t = 0. 

FIG. 4. The time derivative of the index of conformality /conf for the CFL prescription applied on constant-t hypersurfaces, 
for a particle that starts from rest at rp — 17.64(2i\f) [rp = 15(2M)]. The flgure shows that the CFL prescription, in contrast 
to the convective prescription, implies a nonvanishing /conf. 

FIG. 5. Radiated quadrupole energy for a particle falling from rest at ro = 15(2M). The energy actually radiated (labeled 
"True energy" is the horizontal line at E2 = 1.64 x 10~^mo/(2M). The dashed line shows energy computed from convective, 
conformally flat, horizon frozen boundary conditions. The prescribed data was used for a hypersurface on which the particle 
was at position r^, and evolved forward in time. To underscore the role of the Zerilli potential, a scaled plot of the potential is 
shown as a dotted curve. This plot of the potential is shifted vertically so that aligns with the energy curve as r* — » 00. 

FIG. 6. Comparison of the evolved i/> with the new prescription (conformally flat, frozen-horizon, convective) and with BY 
data, on three successive hypersurfaces, for a particle started from rest at ro = 1.5(2M). Even at very early stages the BY 
prescription fails. By contrast, the new prescription continues to give good predictions for most of the trajectory. 

FIG. 7. Comparison of the evolved ?/> with the new prescription (conformally flat, frozen-horizon, convective) and with BY 
data, on three successive hypersurfaces, for a particle started from rest at ro = 15(2M). At very early stages the BY and new 
prescriptions give reasonable accuracy. At somewhat later times BY predictions are greatly in error, while the new prescription 
continues to give good predictions for most of the trajectory outside the peak of the potential. 
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FIG. 8. The waveform as seen by an observer located at 1000(2Af). The solid line is the waveform generated by the infall of 
a particle from rest at ro — 15(2M). The other curves are the result of replacing the evolved data with new data (conformally 
flat, frozen-horizon, convective) on hypersurfaces at the times indicated. If the replacement is made early enough (before the 
particle reaches the maximum of the Zerilli potential) the excitation of quasinormal ringing has the correct magnitude and 
phase, as can be seen in (a). For a replacement at a later time, shown in (6), there is a shift in the quasinormal ringing. 

FIG. 9. The £ = 2 radiated energy as predicted by the extension of the close limit (Ref. [Q) to the particle case. Quadrupolar 
energy calculated from the close limit, applied to the new prescription for data on hypersurfaces at various particle position r^, 
are compared with the true radiated energy. These close limit results are reasonably accurate at sufficiently small separation. 
This should be contrasted with the situation for CFL data, given in Fig. 26 of Paper II, where disagreement was divergent at 
small separation 
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